Weak selection and stability of localized distributions in Ostwald ripening 
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We support and generalize a weak selection rule predicted recently for the self-similar asymptotics 
of the distribution function (DF) in the zero- volume- fraction limit of Ostwald ripening (OR). An 
asymptotic perturbation theory is developed that, when combined with an exact invariance property 
of the system, yields the selection rule, predicts a power-law convergence towards the selected self- 
similar DF and agrees well with our numerical simulations for the interface- and diffusion-controlled 
OR. 
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In a late stage of a first-order phase transition a two- 
phase mixture undergoes coarsening, or Ostwald ripening 
(OR), when the minority phase tends to minimize its 
interfacial energy under condition of a constant volume 
|]-f|. Despite numerous works, OR continues to attract 
attention both in experiment pj and in theory . Our 
main motivation in studying this problem has been an 
attempt to resolve an old selection problem (described 
below) that created much controversy. 

The "classical" formulation of the problem of OR, valid 
in the limit of a negligibly small volume fraction of the 
minority domains, is due to Lifshitz and Slyozov (LS) 
jl],U and Wagner ||]. In this formulation, the dynamics 
of the distribution function (DF) F(R,t) of the domain 
sizes are governed (in scaled variables) by the continuity 
equation 

r)F r) 1 / 1 1 \ 

where R c (t) is the critical radius for expansion/shrinkage 
of an individual drop, while n is determined by the mass 
transfer mechanism. The dynamics are constrained by 
conservation of the total volume of the minority domains: 

/>oo 

/ R 3 F(R, t)dR = Q = const . (2) 
Jo 

Of great interest are possible self-similar intermedi- 
ate asymptotics of this problem and the rule that se- 
lects the relevant asymptotics out of many possibilities. 
Scaling analysis of Eqs. (0) and (||) yields a similarity 
ansatz F(R,t) = $ (Rt~ u ) and R c = (t/af ', where 
(i = 4/(n + 2), v — l/(n + 2) and a — const. Upon 
substitution, one finds a family of self-similar DFs for 
every n > — 1, where each of the DFs is localized on a 
finite interval [0, u m ] of the similarity variable u = Rt~ v . 
The DFs can be parameterized by a, and the interval of 
possible values of a is determined by the requirements 
of the continuity of F{R, t) on the whole interval [0, u m ) 
and normalizability with respect to Eq. (|2|) . For each of 
the solutions, the average domain radius grows in time 
like t 1 /(™+ 2 ), and the number of domains decreases like 
£-3/(n+2)^ ^ut the coefficients in these scaling laws are 
cr-dependent. 



The self-similarity and related scalings were discovered 
by LS |,|| in the case of n = 1 (diffusion-controlled OR) . 
LS arrived at a unique self-similar DF (we will call it 
the limiting solution), and ruled out other possible so- 
lutions. In the first paper JO], the other solutions were 
rejected as non-normalizable with respect to Eq. (|^). 
In the case of n = (interface-controlled OR) this ar- 
gument was repeated by Wagner However, already 
in their second paper [gj on the same subject LS real- 
ized that no problem with normalization arises for ini- 
tially localized DFs (that is, for those with a compact 
support at t = 0). This correction was apparently over- 
looked in the literature {e.g., Ref. 0), until Brown || 
addressed the other solutions and found them numeri- 
cally for n — 1. This created a long-standing contro- 
versy (see, e.g. ||), and the first step towards resolving 
it was made in the case of n = Q. It was noted that 
a DF, initially localized on an interval [0, R m {t = 0)], 
always remains localized on a (time-dependent) interval 
[0, Rm(t)]. Furthermore, if F(R, t = 0) is describable by 
a power law Aq [R m (t = 0) — R] x in the close vicinity of 
R = R m {t = 0), then for any t > the leading term in 
the expansion of F(R, t) in the vicinity of R = R m {t) has 
the form A(t) \R m (t) — R] x . Invariance of the exponent 
A under the dynamics (Q) and (^) implies a (weak) selec- 
tion rule for the "correct" self-similar DF, as there is a 
one-to-one correspondence between < A < +oo and the 
parameter a |6| . [The limiting solution corresponds to an 
extended (non-compact) initial condition or, formally, to 
A — > +oo.] More precisely, if a self-similar asymptotics is 
ever reached, it must be the one selected by A. However, 
no attempts have been made to solve (even numerically) 
the full time-dependent problem with a localized initial 
DF. Furthermore, no stability/convergence analysis for 
the localized DFs has been performed, so the selection 
rule proposed in Q has remained unconfirmed. 

This Letter supports the selection rule along three di- 
rections. The first one is to generalize the selection rule 
for any n > — 1. The second is to prove the stability 
of and analyze the convergence towards the selected self- 
similar DF. The third is to verify our theoretical predic- 
tions numerically. 
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A meaningful formulation of the stability problem re- 
quires some care. Indeed, each member of the family of 
self-similar solutions for the DF, except the limiting so- 
lution, is formally unstable with respect to addition of an 
(infinite) tail. In this case it is the limiting solution that 
will finally develop |jl|]|J. However, such a perturbation 
is not always possible. In addition, the results of |gJ im- 
ply that each member of the family, except the limiting 
solution, is formally unstable with respect to a localized 
perturbation that either has a larger R m than the "un- 
perturbed" DF, or the same R m and a smaller exponent 
A. In each of these cases another self-similar solution 
from the same family finally develops (as we see in our 
numerical simulations), and this situation can hardly be 
regarded as instability. A meaningful formulation of the 
stability problem should therefore deal with initial per- 
turbations localized on the same interval of R as the "un- 
perturbed" DF, and characterized by the same exponent 
in the close vicinity of R m (t = 0). 

We will develop an asymptotic linear theory that, com- 
bined with an (exact) invariance property of the model, 
will enable us to prove, analytically, the stability of each 
of the self-similar DFs. This result and our numerical 
simulations will confirm the weak selection rule 0. We 
will analyze the late-time convergence of an initially lo- 
calized DF towards the selected self-similar DF and find a 
power-law decay in time for the corresponding (non-self- 
similar) perturbation. This decay is much faster than the 
logarithmic decay found for the limiting solution KB. 
We will see that not only the selected self-similar DF, 
but also the decay exponent are determined solely by the 
analytical properties of F(R, t = 0) in the close vicinity 
of R m (t — 0). Our theoretical predictions show a very 
good agreement with simulations. 

We will start with the asymptotic theory. Solving the 
problem analytically is made possible by a change of 
variables that employs the compactness of the support 
[0, R m {t)} of the DF. Introduce a scaled drop radius and 
a new time: 



x(R,t) 



R 



Rm (t) 



and r 



dt' 



R m +2 (t>) ' 



(3) 



and a scaled DF 



G(x,T)=Ri[t(T)}F[R(x,r),t(r)}. (4) 

In the new variables Eqs. ([!]) and (|J) can be rewritten as 

(dG/dr) + [v (x~ n -x)+x- x-"- 1 ] (dG/dx) 

+ [(n + l)x- n - 2 - nvx- 71 - 1 - A(v - 1)] G = (5) 

l 

and J G(x,t)x 3 dx = Q, respectively, where v(t) = 
o 

i? m [t(r)]/ R c [t(r)]. The function G(x,t) is nonzero on 
the interval < x < 1 and zero elsewhere. 



We will see in a moment that a self-similar solution for 
F(R, t) corresponds to a steady-state solution for G{x, r). 
Therefore, we are looking for the solution in the following 
form: 



G{x,t) = $o(x) + $i{x) e qT + 
v(t) = v + V\ e qT + 



(6) 



where q is a (sought for) complex number. Both <&o(x), 
and <f>±(x) are localized on the interval [0, 1]. The pertur- 
bation must not change the normalization condition (|^) 

and the asymptotics of the unperturbed solution near the 
l 

point x = 1, that is, / x 3 <&i(x) dx — and <&i = O($ ) 
o 

at x — > 1. 

A family of steady-state solutions &o(x) (parameter- 
ized by vq) is obtained from the zero-order equation 



[vq (x n — x) + X — X 



(d$ /dx) 



+ [(n+l)x 



-n-2 



UVqX 



-n-1 



4(«o-l)] $o = 0. (7) 



Integration of this equation in elementary functions is 
possible for n = — 1, 0, 1 and 2 (that is, for most cases of 
physical interest jn|). For example, for n = one has 

$ (x) = C Q x (1 - x) a {x 2 - x) v , < x < 1 , (8) 

where 



4^o-5 
2-Uo 



v o - 5 
2-«o ' 



■r-2 



Vq-1 



(9) 



while Cq is determined from the condition J x 3 $o (%) dx 

o 

= Q. This family of solutions is defined for 5/4 < vq < 2. 
It corresponds to the family of self-similar solutions for 
F(R, t) obtained in Ref. § . 
For n = 1 one obtains 

$o = Cqx 2 {1 - x) a (x - x-y 2 -^(x+ - x)" 11 -' 12 . (10) 

Here 



5t>o — 6 
3-2vo 



7i 



12 - 7v 
6-4u 



72 



3v 



(6 - 4v )s ' 

(11) 



x± = (-l±s)/2, s = [(vq + 3)/v -1)} 1 / 2 andO < x < 1. 
This family is defined for 6/5 < vq < 3/2. 

For any n, we will need to know the behavior of $o (%) 
and $i (x) in the close vicinity of x = 1. A simple analysis 
of Eq. (0) yields $ = flo(C) C , where C = 1 - x, #o(C) 
is an analytic function on the interval [0,1], i?o(0) ^ 0, 
and 



(n + 4) vq — n — 5 
n + 2 - (n + 1) v 



(12) 
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The solution for &q(x) exists if < a < oo, that is 
(n + 5)/(n + 4) < v < (n + 2)/(n + 1). This interval 
of permitted values of vq is non-empty for any n > —1. 
[The case of n = — 1 is the simplest: Hq(() — const.] 
Now we go to the first order in Eq. (||): 

[v {x~ n -x)+x- x-' 1 - 1 ] {d^i/dx) 
+ [q-(n+ l)x- n - 2 - nvox-™- 1 - 4(u - 1)] $i 

= vi [{x - x~ n ) {d&o/dx) + (4 + k-"- 1 )^] . (13) 

For a given v\, this linear equation can be solved in 
quadratures ]Tl| ]. We will need only the leading asymp- 
totics of the solution in the close vicinity of x = 1 , so we 
write down the solution as 



$i = Cxi(C) +C /3 x 2 (C), 



(14) 



where (3 = a — q[n + 2 — (n + l)i> ] -1 , Xi and X2 are 
analytic functions on the interval [0, 1] and Xi,2 (0) 7^ 0. 
The solution exists if Re (3 > which implies Re q < 
(71 + 4)«o — n — 5. One can check a posteriori that this 
inequality holds. 

Eq. ( |l4| ) will be used later. At this stage we notice 
that the still undetermined "eigenvalue" q must be se- 
lected by the initial condition. To make this selection 
possible, we should exploit an (exact) invariance prop- 
erty of Eq. (Q). Consider the initial value problem 
dR/dt — V(R,t),R(0) — Ro that describes the char- 
acteristics of Eq. (Q). If the solution of this problem, 
R(t;Ro), is known, the solution of Eq. ([TJ) can be writ- 
ten in the following form: 



F(R,t) = F [R {R,t)]dR (R,t)/dR, 



(15) 



where Ro(R, t) is the function inverse to R(t, Rq) with re- 
spect to the argument Rq. Now, R{t, Ro) is an analytic 
and monotonic function of Rq. Therefore, the inverse 
function Rq (R, t) is also an analytic function of R, and 
so F(R, t) preserves its analytic form along the character- 
istics R = i?(t;i?o), including the "edge" characteristics 
R m (t). 

We assume a power-law behavior of F(R 7 t = 0) in the 
close vicinity of R = R m (0). More precisely, we assume 
that 



(16) 



where £ = i? m (0) — R > 0. Here Ai and A2 > Ai are ar- 
bitrary positive numbers, such that A2 — Ai 7^ 1, 2, . . . , 
and gi(£) and (72(C) are analytic at £ = such that 
31,2(0) 7^ 0. In view of of the analyticity property men- 
tioned above, Eq. (15) can be rewritten as 



F(R,t) = (£') Al Mr,i) + (0 A2 ^(£',i) , (17) 

where £' = R m (t) - R > 0, while fa(£',t) and fa(£',t) 
are analytic functions of £' at £' = 0, and fa, 2(0, t) 7^ 0. 



Under the transformation (||), ([IJ) the variables i? and 
F are multiplied by some quantities independent of R. 
Therefore, we can rewrite Eqs. ( fHf ) and (|l7|) in the new 
variables x and r as follows. The initial DF is now 



Go(x)=Ci A ffi(0 + C A2 3 2 (0, C>0, 



(18) 



where g[(C) and g' 2 (() are analytic at ( — 0, 2(0) 7^ 0> 
and we remind that £ = 1 — x. Correspondingly, the 
time-dependent DF G(x, r) can be written as 

G(x,T)=( x iti 1 ((,T)+( x *ti 2 ((,T), C>0, (19) 

where h\{Q, r) and /j 2 (Cj t ) are analytic functions of £ in 
C = 0, and fa. 2 (0,r) ^ 0. 

The exponents Ai and A 2 , prescribed by the initial con- 
ditions, remain invariant. Hence, the long-time asymp- 
totics of Eq. (|l9|) should coincide with that given by 
Eqs. (||) and (|14|). A direct comparison yields a = X\ 
and (3 = A 2 . The first equality is nothing but a (weak) se- 
lection rule for the self-similar solution, and the selected 
value of vq is 



^0 



(n + 2)Ai +n + 5 
(rc+l)Ai + n + 4 

The second equality determines q: 

3(A 2 -A0 



q = 



(n + l)Ai +n + 4 



(20) 



(21) 



One can see that —q is real and positive which means 
stability. 

Going back to the "physical" variables R and t is easy. 
Indeed, evaluating R m (t) on the self-similar solution, we 
obtain R m (t) = [(n + 2) (u - 1) t] 1/{n+2) . Then, using 
Eq. (|J), we see that e 9T = t~ r , a power-law decay in the 
"physical" time. Here 



3(A 2 -A!) 
(n + 2)(Ai + : 



> 0. 



If we limit ourselves to an important particular case 
of a single "non-trivial" exponent in the initial DF, 
Gq(x) = ( x g((), [where g(() is analytic and 5(0) 7^ 0], 
then G(x,t) — £ a /i(£,t), where faCi 7 ") is an analytic 
functions of ( and /i(0,r) 7^ 0. Now, using Eqs. (||) and 
(|l4|), we obtain 



q 



(n + 1)A + n + 4 



(22) 



unless the linear term in the Taylor series of g(() at ( = 
vanishes |fl2| . This yields the power exponent 



(23) 



(n + 2)(A + l) 



In the limit of A — > 00 , we obtain T — > 0. Clearly, it 
corresponds to the logarithmically slow decay obtained 
for the limiting solution |l|,||. 
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Therefore, both the self-similar DF, and the power-law 
decay rate of a small perturbation around it are uniquely 
determined by the asymptotics of the initial DF in the 
close vicinity of the maximum domain size, R = R m (0). 

We verified the theory (in the cases n — and 1) by 
performing extensive numerical simulations with Eq. (jl|) 
and an explicit equation for R c that follows from Eqs. 
(|l|) and (Q). As the dynamics is extremely sensitive to 
small changes in the vicinity of R = R m (t), we needed an 
algorithm that preserved the compactness of the DF and 
kept a high accuracy near the edge point R = R m (t). A 
simple and efficient Lagrangian algorithm was developed 
fl3f that satisfied these requirements. Typical simula- 
tion results for the interface-controlled OR, n = 0, are 
presented in Figs. 1 and 2. Figure 1 shows convergence 
of an initially localized DF, F(i?,0) = R (5 - R) x with 
A = 1, towards the selected self-similar DF (||), for which 
Eq. ( po|) predicts v — 7/5. The inset shows convergence 
of v(t) towards vq = 7/5. The convergence exponent 
Teiep = 0.76 found numerically agrees very well with our 
theoretical prediction Y t h = 0.75. Figure 2 shows the 
convergence exponents T exp found numerically for dif- 
ferent A. A good agreement with the theoretical curve 
Tt/i = 3/[2(A+l)] is seen. We also observed a good agree- 
ment between the theory and simulations in the case of 
the diffusion-controlled OR, n = 1. 
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FIG. lP-0 Convince oPa4i initial-^ localized? DF wM> 
A = 1 towards the selected self-similar DF (^) with vo = 7/5 
(solid line). Numerical solutions are shown by dotted lines 
at time moments t = 20 (a), 100 (b), 500 (c) and 1000. The 
inset shows the convergence of v(t) towards Vo = 7/5. 

We have demonstrated that only a weak selection is 
possible in the "classical" model of OR. To get a strong 
selection rule, one obviously must go beyond the "classi- 
cal" model. 
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FIG. 2^ fic \Gonvergenefivexponen1is!^predictedlalnalytically 
(line) and found numerically (squares) for different A. 
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